Magneto-optical study of magnetic flux penetration into a current-carrying 
high-temperature superconductor strip 
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Magnetic flux distribution across a high-temperature superconductor strip is measured using 
magneto-optical imaging at 15 K. Both the current-carrying state and the remanent state after 
transport current are studied up to the currents 0.97/c where Ic is the critical current. To avoid 
overheating of the sample current pulses with duration 50 ms were employed. The results are 
compared with predictions of the Bean model for the thin strip geometry. In the current-carrying 
state, reasonable agreement is found. However, there is a systematic deviation - the flux penetration 
is deeper than theoretically predicted. A much better agreement is achieved by accounting for flux 
creep as shown by our computer simulations. In the remanent state, the Bean model fails to 
explain the experimental results. The results for the currents / < 0.7/c can be understood within 
the framework of our flux creep simulations. However, after the currents / > Q.77c the total flux 
trapped in a strip is substantially less than predicted by the simulations. Furthermore, it decreases 
with increasing current. Excessive dissipation of power in the annihilation zone formed in the 
remanent state is believed to be the source of this unexpected behavior. 

PACS numbers: 74.76.Bz, 74.60. Ge, 74.60. Jg,78.2QLs 
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I. INTRODUCTION 

Spatially resolved studies of magnetic flux penetration 
into high-temperature superconductor (HTSC) films are 
extensively performed during the last few years. Modern 
experimental techniques, in particular Hall micro-probe 
measurements and magneto-optical (MO) techniques, al- 
low local magnetic field distributions in various HTSC 
structures to be investigated with rather high spatial res- 
olution. As a result, a quantitative comparison between 
experimental flux density profiles and theoretical predic- 
tions has become possible. Most of the comparisons are 
done within the framework of the critical state model 
(CSM)EJ. Based on this model theoretical calculations of 
the field distributions for many practical geometries are 
carried out. In particular, the flux profiles in an infinite 
thin strip placed in a perpendicular magnetic field, or 
carrying a transport current are calculated in Refs. |^-^. 
Most experimental studies of fiux penetration, see e. g. 
Refs. §-| are focused on the behavior of samples placed 
in an external magnetic field, as well as on the remanent 
state after the field is switched-off. The results appear in 
a good agreement with experiment. 

Meanwhile, we are aware of only a few papers devoted 
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to experipiental investigation of the self-field of transport 
currentaJ^Ej. Unfortunately, these investigations do not 
allow ai^iiaple comparison to the theory. Indeed 
of thenJj'0~E3 present results for samples of pit-tif 
plicated geometry, e. g., for tapes. 0therg|l3O report 
results for currents much less than the critical current, 

Ic 

The aim of this work is to study by MO imaging the 
flux penetration into a strip with transport current. It 
seems most interesting to have the current / as close as 
possible to the critical one. To reach this aim one needs 
narrow strips where the critical current is not too large to 
be carried by contacts without their destruction. On the 
other hand, the wider the strip the better the relative spa- 
tial resolution of MO imaging. By optimizing both the 
strip's width and other experimental conditions we man- 
aged to obtain fiux profiles both in the current-carrying 
and in the remanent state with a resolution sufficient for 
comparison with theory. 

In Sec. H the samples and the experimental technique 
are described. The results for flux proflles and recon- 
structed current distributions are compared to the CSM 
in Sec. III. It is shown that the deviations are fairly 
small in the current-carrying state. However, they are 
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pronounced in the remanent state after applied current. 
The main deviation is a deeper flux penetration inside the 
strip compared to predictions of the CSM. To understand 
the source of the deviation we have carried out numer- 
ical simulations of the field and current profiles taking 
into account flux creep. The results of these simulations 
are discussed in Sec. IV. Creep appears to explain all 
the experimental results for the current-carrying state, 
as well as the results for the remanent state after rela- 
tively small currents, / ^ 0.7 /c. In the remanent state 
after large currents the experimental profiles could not 
be explained either by CSM or by flux creep. It seems 
that thermal effects are responsible for this behavior. 



II. EXPERIMENTAL 
A. Sample preparation 

Films ot^Ba2Cu307_5 were grown by dc magnetron 
sputteringl3 on LaAlOa substrates. X-ray and Raman 
spectroscopy analysis confirmed that the films were c- 
axis oriented and of a high structural perfection. Several 
bridges were formed from each film by a standard lithog- 
raphy procedure. Their dimensions are 500 x 110 x 0.2 
/im^. To minimize temperature increase caused by Joule 
heating the contact pads were made as wide as possible. 
They are displaced to the side of the structure allowing 
the MO indicator film to be placed as close to the bridge 
as possible. To provide low contact resistance they were 
covered with a Ag layer, and Au wires of 50 fim diame- 
ter were attached by thermal compression. The bound- 
ary resistance of Ag/YBa2Cu307_5 interface was as low 
as 10"'^ — 10"^ ft crn^, while the area of one contact pad 
was Ri 0.25 cm^. The bottom of 500 micron substrate 
was held at a constant teiaperature. Since the thermal 
conductivity of LaAlOs islla 0.1 W/ (cm K), the temper- 
ature rise at the interface induced by Joule heating due 
to currents up to 6 A is always less than 0.1 K. Thus, re- 
sulting heating of the YBa2Cu307_5 bridge situated 500 
micron away from the contact pads is negligible. 

An initial selection of the structures having the 
smoothest surface, i.e., height of over-growth less than 
3 /im, was made using SEM. Next, bridges with pro- 
nounced weak links, Tc-inhomogeneities and other de- 
fects which reduce the total critical current Ic were elimi- 
nated by means of low-temperature SEM,Ej in addition to 
current-voltage measurements made by a standard four- 
probe scheme. As a result, the bridges used for the final 
investigations had a critical current density, jc, larger 
than 10^ A/cm^ at 77 K. 

B. Magneto-optical imaging 

Our fiux visualization system is based on the Faraday 
rotation of a polarized light beam illuminating an MO- 



active indicator film that we place directly on top of the 
sample's surface. The rotation angle increases with the 
magnitude of the local magnetic field perpendicular to 
the HTSC film. By using crossed polarizers in an optical 
microscope one can directly visualize and quantify the 
field distribution across the sample area. As Faraday- 
active indicator we use a Bi-jJoped yttrium iron garnet 
film with in-plane anisotropycj. The indicator film was 
deposited to a thickness of 5 /im by liquid phase epitaxy 
on a gadolinium gallium garnet substrate. Finally, a thin 
layer of aluminum was evaporated onto the film in order 
to reflect the incident light and thus providing a double 
Faraday rotation of the light beam. The images were 
recorded with an eight-bit Kodak DCS 420 CCD camera 
and transferred to a computer for processing. The con- 
version the gray level of the image into magnetic field val- 
ues is based on a careful calibration of Bi:YIG indicator 
response to a range of controlled perpendicular magnetic 
field as seen by the CCD camera through the microscope 
(see also Ref. ^). After each series of measurement, the 
temperature was increased above Tc and the in-situ cali- 
bration of the indicator film was carried out. As a result, 
possible errors caused by inhomogeneities of both indica- 
tor film and light intensity were excluded. 

To avoid overheating of the HTSC bridges as the cur- 
rent approaches Ic a specific experimental procedure was 
developed. Short current pulses were applied and syn- 
chronized with the camera recording as shown in Fig. n^. 
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FIG. 1. Temporal profile of the pulsed transport current. 
The exact time intervals of the MO image recording in the 
current-carrying state and the subsequent remanent state are 
indicated. 

The transport current with a rise time of 10 ms was ap- 
plied 20-30 ms before an image was recorded. The expo- 
sure time was 35 ms, after which the current was ramped 
to zero in 10 ms. After waiting another 20-30 seconds the 
remanent field distribution was measured. 
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III. RESULTS AND COMPARISON WITH CSM 

Throughout the paper we use the following notations, 
see Fig. g. The a:-axis is directed across the bridge, the 
edges being located at a; = ±w. The y-axis points along 
the bridge, and the z-axis is normal to the film plane. The 
z-component of the flux density is denoted by B. The 
profiles were always measured for a fixed y in the 

central part of the strip minimizing the effect of the stray 
field of the contact pads. Although the MO measure- 
ments give the distribution of we could in the sim- 
ple geometry under consideration always determine the 
sign of B by inspection of the images. We let J{x) denote 
the sheet current density defined as J{x) = J j(x,z)dz, 
where j{x,z) is the current density. For brevity, we will 
often use current density also for J{x). As the bridge 
thickness is much less than its widLh, the theoretical re- 
sults for the thin strip geometryoQ are used below. In 
our experiments, the bridge thickness is also of the order 
of the penetration depth A, hence its magnetic proper- 
ties are fully characterized by the two-dimensional flux 
distribution at the surface. 



A. Current-carrying bridge 




FIG. 2. (a) Sketch of the studied YBa2Cu307_i bridge 
structure, (b) MO image of flux distribution with the bridge 
in the current-carrying state with / = 4.16 A. (c) MO image 
of flux distribution with the bridge in the remanent state after 
the applied-current state shown in (b) . Two inner bright lines 
near the sample edges represent the original trapped flux, the 
outer bright lines represent penetration of the oppositely di- 
rected return field, while dark lines between them correspond 
to regions of vortex-antivortex annihilation. Strong contrast 
enhancement was applied to the image in (c). 



Flux density distributions for a strip carrying a 
transport current / were measured for currents up to 
7=5.72 A. The MO image for the current / = 4.16 A 
is shown in Fig. ^(b). Three profiles of the flux density 
taken across the strip are shown in Fig. ||. The 
profiles have maxima near the strip edges and a minimum 
in between. Actually, the left and the right parts of the 
profiles correspond to induction of opposite sign. As the 
current increases, (a)-(c), the magnetic flux penetrates 
deeper and the flux free Meissner region in the center of 
the bridge decreases in size. 

The experimental MO data are interpreted in the 
framework of the CSM. As discussed in Ref. H, one should 
account for the finite distance between the sample and 
the MO indicator film. The perpendicular magnetic field, 
B, at the height h above the center of the bridge can b| 
calculated from the current density distribution J{x) at 
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w + {x' — x)' 



J{x')dx' +Ba. (1) 



Here Ba is the external magnetic induction. The cur- 
rent density distribution in a strip carrying a tr|ansport 
current / can be written for the Bean model asau 



■M. 



— arctan 

TT 













-X2 1 



\x\ < a 



(2) 



a < |a;| < w 



3 



where a = wy/1 — (J/Ic)'^, and = 2w is the critical 
current. As Eqs. (0) and (Q) yield a symmetric |i3(a;)|- 
profile we find it necessary to account also for the stray 
field of the contact pads in order to reproduce the slight 
asymmetry in the |i?(a;)| data. The current in the pads 
produces near the central part of the bridge a magnetic 
field which acts as an additional external field varying 
slowly in space. This allows us to employ the results 
of the CSM for the case of a transport cucrent super- 
imposed by SLyMeak external magnetic field.Eil According 
to this theoryOu, the current density distribution Eq. (g) 
is modified to yield J{x) = Jc at x within the intervals 
{~w,p~a), (p + a, w) and 



J{x) 



if , (x — p)(w — p) — a 
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a{w + x) 
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for p — a<x<p + a. Here 



cosh{Ba/Bc) 
w— tanh — - 

L \ Br 



I 

z 

Br = 



100 

50 




Profiles of |i?(a;)| calculated from Eqs. (^) and (P are 
shown in Fig. ^ by the solid lines. Here Ic and h are 
parameters determined by fitting the CSM profiles to ex- 
perimental data. The values Ic = 5.9 A, h = 8 fim yield 
the best fit for the profiles measured at all currents /. 
Note that Ic has a physical meaning of the critical cur- 
rent for CSM model, at which the magnetic field fully 
penetrates the bridge. 

One can notice in such experiments, and e. g. in 
Fig. |^(b) , that the experimental penetration of the mag- 
netic field is deeper than the CSM prediction. This char- 
acteristic deviation can be accounted for by introducing 
flux creep, as will be discussed in detail in the next sec- 
tion. 

With h being a known parameter one can also directly 
determine the current distributions across the strip, J{x), 
from the experimental _B(x)-profiles. This can be done 
on a model-independent basis according to an inversion 
scheme developed in Ref. ||. Since the accuracy of the in- 
version is limited by the distance /i, we chose a step size 
along the a;-axis close to A = h/2. Specifying the coor- 
dinates in units of A, i.e., x = nA, x' = n' A, h — dA, 
and applying a Hannig window filtering function one 
obtainsEl: 
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FIG. 3. Profiles of the perpendicular magnetic field pro- 
duced at a distance h = 8 nm above the bridge carrying trans- 
port currents. Symbols show the experimental data (only 1/5 
of all measured points are shown) , and the solid lines represent 
CSM calculations from the Eqs.(|) and @. I/h = 0.13(a), 
0.53(b), 0.97(c), where 7c = 5.9 A. 



The current density profiles calculated from experi- 
mental -B(n)-data using this formula are shown in Fig. 0. 
The figure also shows the results obtained from Eq. (|^) 
using the value of Ic determined by the fitting of |_B(a;)|- 
profiles. As seen from Fig. ^, the experimental profiles 
trace all qualitative features of the theoretical curves. 
In particular, the minimum in the current density near 
the bridge center predicted by the CSM can clearly be 
distinguished in all experimental curves. Moreover, the 
position of the minimum is found to be slightly shifted 
towards negative X, Oil behavior in full agreement with 
the theory when the effect of contact pad stray fields is 
accounted for. 

As first pointed out in Ref. |[ the MO YIG indicators 
with in-plane anisotropy respond not only to Bz but also 
to the component Bx parallel to the film. In the data pre- 
sented we always made the proper corrections according 
to the method suggested in Ref. ||, i.e., Bx{x) component 
was calculated from the current profiles by Biot-Savart 
law and then used to re- normalize the profile Bz{x). This 
procedure is converging very rapidly. 

To check self-consistency of our inversion calculations 
we integrated the current density over x. Indeed, The 
total current was always equal to the transport current 
passed through the bridge within an accuracy of 5%. 
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FIG. 4. Distributions of current density in the cur- 
rent-carrying strip. Symbols show the distributions inferred 
from experimental data using Eq. (^). The predictions of the 
GSM, Eq. (|) are shown as full lines. I/h = 0.13(a), 0.53(b), 
0.97(c), where h = 5.9 A. 

Some deviations between experimental and theoretical 
current profile slopes are seen near the sample edges 
where the theoretical profiles are discontinuous. The ob- 
served smearing of the current profiles reconstructed from 
the induction data can be ascribed mainly to discreteness 



of the points used in the inverse calculation, and possibly 
a slight structural imperfections near the very edge of the 
sample. 



B. Bridge in a remanent state after transport 
current 

When the transport current / is switched-ofF, the mag- 
netic flux in the inner part of the strip remains trapped. 
The return field of this trapped flux will re-magnetize the 
edge regions of the strip and the flux of opposite sign pen- 
etrates an outer rim. As a result, the measured remanent 
distribution |-B(a::)| should display two peaks in each half 
of the bridge; one representing the maximum trapped 
flux, and another near the edge indicating the maxima 
in the reverse flux. In this oscillatory flux profile typical 
values of \Bz\ are significantly lower than in the current- 
carrying state. Thus, MO studies become substantially 
more difficult in the remanent state, and flux profiles with 
reasonable signal-to-noise ratio were recorded only after 
relatively large transport currents, / > 0.6/c = 3.68 A. 

The MO image of flux density distribution in the rema- 
nent state after switching off a transport current of 4.16 
A is shown in Fig. ^c). Figure ^ shows the flux density 
profile taken across the strip in this remanent state. It 
can be seen that there are large maxima of trapped flux 
in the center of the bridge. Also, the weaker maxima of 
reverse flux are visible near the edges. In the figure we 
have removed the experimental points corresponding to 
the regions where the MO image is governed by disturb- 
ing zig-zag domain pattern in the indicator film. This 
occurs where changes sign. 

The current distributi(Mj-.derived from the Bean model 
for the remanent state iscia 
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-a < X < a, 



a <\x\ <b 



(5) 



where a = w^l - {I/IcY, b = w^/l - {I/^hf and / is 
the maximal current. At 6 < \x\ < w the current den- 
sity is equal to —Jc- The small contribution from the 
field generated by the contact remanent currents is here 
neglected. 

The magnetic field profiles calculated using Eqs. (|^) 
and (|l|) are shown in Fig. || (solid line) together with the 
experimental data. We used the same values for Ic and h 
as obtained by fitting results for current-carrying state. 



Evidently, there is here a significant deviation between 
our data and the CSM description. The main devia- 
tion is that the trapped flux maxima in the experimental 
curve are shifted towards the center. Again, this can be 
shown to be an effect of flux creep as discussed in the 
next section. Another deviation seen in Fig. |^ is that 
the maxima near the edges are hardly visible experimen- 
tally. This is most likely due to the nonlinear response 
of the optical detection system which leads to reduced 
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sensitivity at low Faraday rotation angles, i.e., at small 
induction values. Another source of smearing is that the 
width of the maxima near the edges is comparable to the 
thickness of MO indicator film. 

The experimental flux profiles in the remanent state 
after large currents / > 4.16 show even more intrigu- 
ing behavior. Fig. ^ shows six remanent flux profiles 
across the right half of the bridge after applying cur- 
rents ranging from 4.42 A to 5.72 A. The observed de- 
crease of the trapped flux with increasing transport cur- 
rent seems highly unexpected. Such a behavior definitely 
contradicts the CSM. To make a quantitative compar- 
ison we have integrated the flux trapped in the band 
— 0.4?ii < X < Q.2w , where the MO data are reliable. 
The dependence of the total flux on the current is shown 
in Fig. ^. At small currents, / < 4.16 A, the behavior is 
normal as the amount of trapped flux increases with the 
current. At larger currents, however, the trapped flux 
levels out and even starts to decrease as I approaches 7c- 
The full line in the same figure indicates the predictions of 
the CSM, which clearly shows a different behavior. Nei- 
ther is the non-monotonous behavior seen experimentally 
explained by flux creep. 

Let us sum up the comparison with the CSM for the 
current-carrying and remanent states. The general fea- 
tures of flux distribution in a strip with transport current 
predicted by the CSM are well confirmed by the exper- 
iments. It should be noted that such a good agreement 
was achieved only by taking into account a finite distance 
between the HTSC film and the MO indicator, the 




FIG. 5. Remanent profiles of the absolute value of per- 
pendicular magnetic field at the distance h = 8 fim above the 
bridge after applying transport current / — 4.16 A. Dots - 
experiment, solid line - CSM, Eqs. (|), (Q) with Ic = 5.9 A, 
dashed line - flux creep simulations. A large amount of flux 
trapped near the bridge center is in strong contradiction with 
CSM and can be explained by flux creep. Experimental points 
in the regions where the parallel component Bxix) changes 
sign and the response of indicator film is suppressed have 
been removed from the plot. 




FIG. 6. Experimental profiles of the absolute value of per- 
pendicular magnetic field in the remanent state after different 
transport currents I: (1) - 4.42 A, (2) - 4.68 A, (3) - 4.94 A, 
(4) - 5.20 A, (5) - 5.46 A, (6) - 5.72 A. Only the right haff of 
the bridge is shown. Note the unusual /-dependence of the 
trapped fiux: the larger the current the smaller the density 
of trapped fiux in the bridge. 
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FIG. 7. Flux trapped within the band —QAw < x < 0.2w 
in the remanent state after transport current / as a function 
of I. The band corresponds to reliable MO data shown in 
Fig. ^ the fiux in normalized to the band width. Dots - ex- 
perimental results. Solid and dashed lines - CSM and flux 
creep simulations respectively. For 7 < 4.16 A the experi- 
ment is in a qualitative agreement with flux creep, however 
for larger currents a striking steady decrease in trapped flux 
is observed. 

additional field generated by contact currents, and the 
influence of parallel field component on the properties 
of YIG indicator. However, a general trend of a deeper 
flux penetration compared to the CSM prediction can 
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be traced. In the remanent state the deviations from 
the CSM are more pronounced, and experimental flux 
profiles after currents / > 0.71c have qualitatively dif- 
ferent shape. In particular, the flux trapped inside the 
bridge can be several times greater than predicted and 
it depends on the previously applied current in a non- 
monotonous way. 

Comparing the experimental results with the predic- 
tions of CSM we employed the expressions (^) and (j^) 
based on the Bean model. Thus, we assumed the critical 
current density Jc to be _B-independent. We have checked 
validity of this assumption studying penetration of an ap- 
plied magnetic field into a bridge of the same geometry 
fabricated on the same fllm. We could trace rather weak 
dependence of Jc versus B only above B ^ 200 G. For 
smaller B all the results, including the results for the re- 
manent state after external field, were compatible with 
the CSM predictions, as in Ref. |5[ Consequently, there 
is no reason to attribute the deviations observed in the 
transport current regime to a i?-dependence of Jc- 



IV. FLUX CREEP 

The basic assumption of the CSM is that in the re- 
gions where the local current density J is less than the 
critical one Jc, the flux lines do not move. This assump- 
tion does not hold for any finite temperature because of 
thermally-activated flux motion or flux creep. As a re- 
sult, a small amount of magnetic flux will penetrate into 
the regions with J < Jc- We suggest that this excess flux 
penetration is responsible for the smoothening of the ex- 
perimental proflles B{x) relative to the CSM proflles. 

To verify this suggestion we have carried out computer 
simulations of flux penetration into a strip in the flux 
creep regime. Simulations of this type have previously 
proven to be very powerful in the analysis of flux behav- 



ior. In particu 
tization curves 
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ix^hey have been used to anal 
3Ej, magnetic relaxation dataE 
susceptibility for various geometriea^J. Flux creep simu- 
lations have also been used to explain features of the flux 
penetpation into thin HTSC samples in applied magnetic 
fielcB'cSl. For the first time we present here creep simula- 
tions to analyze flux dynamics in a strip with transport 
current and the subsequent remanent state. 



A. Model 

Motion of magnetic flux is governed by the equation 

(6) 



where v is the vortex velocity. The velocity v is assumed 
to be dependent on the current density J and the tem- 
perature T as V = t^o exp[— f/(J)/fcT] where U{J) is the 



current-dependent activation energy due to vortex pin- 
ning. Its dependence on the current density is extensively 
discussed in the literature, see for a review Refs. [27||2^ . 

A conventional approach to the flux creep is based 
upon the linear (Anderson-Kim) relation U{J) = C/c(l — 
J / Jcp) for the pinning energy. Here we use notation Jcp 
for the depinning critical current density which is dif- 
ferent from the CSM critical current density, Jc- The 
Anderson-Kim approximation is fairly good at Jcp — J <C 
Jcp- However, due to low pinning energies in HTSCs, 
the time dependence of the current density J appears 
very pronounced, and during the time of experiment J 
can reach the values well below Jcp- As a result, to de- 
scribe experimental data one should take inta-jaccount 
a nonlinear character of the U{J) dependencd2Zlc3. The 
usual way is to express the pinning energy as a power law 
function of the current density, U{J) oc J~^. Such a de- 
pendence follows from several theoretical models based 
on the concept of collective creep, the values of /i be- 
ing dependent on vortex density, the current, the pinning 
strength, and dimensionalitycZI. Since we are interested in 
the region of low flelds, the value ^ = 1/7 seems most ap- 
propriate, as it applies to the case of a single- vortex creep 
at high current density and low temperaturecj. We use 
the expression U{J) = Uc [{Jcp/ J)^ — 1] for the activa- 
tion energy. In such a normalization U{Jcp) = and Jcp 
retains the meaning of a depinning critical current. The 
current-independent factor in the vortex velocity equals 
vq eKp(Uc/kT) where vq is a quantity of the order of the 
flux-flow velocity. 

The numerical integration of equation was carried 
out by a single step method, similar to the one reported 
in Ref. |2^. Having the field distribution B{x,t) at time 
t, we calculate the corresponding current density distri- 
bution asQ 



J(x,t) = 



B{x' ,t) - Ba{t) w'^-x' 



X — X' 



-dx' 



m 



(7) 



Here I{t) and Ba{t) are the time-dependent total trans- 
port current and applied field, respectively. The expres- 
sion follows from the Maxwell law for the thin film ge- 
ometry. Then the quantity dB{x, t)/dt is calculated from 
Eq. (^). The new distribution of magnetic field is calcu- 
lated as B{x,t) + SB{x,t), where SB{x,t) = 8t-{dB/dt). 
Here a time increment 5t is chosen so that SB(x,t) < 
<5Smax = O-OOOlfXoJcp for any x- Then we come to the 
next step and so on. 

There are two independent parameters in the model 
- the product vq exp{Uc/kT) and the ratio J^pUd [kT)- 
Unfortunately, it seems very difficult to estimate these 
parameters from the literature because of a very large 
scatter in the published values of the pinning energy Uc 
for YBa2Cu307_5. We determine the ratio Jl^pUc/ {kT) 
by fitting the experimental flux profiles and estimate the 
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quantity vq exp(Uc/kT) from a separate experiment on 
magnetic relaxation. 



o 




T 1 1 1 r 
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The experimental quantity In (— (?ln J5/9t) is plotted as 

a function of _B ^ in Fig. ^ The experimental points 
show a nearly linear dependence, and the line repre- 
senting the best fit is also shown. From the fit wc ob- 
tain Wo exp{Uc/kT) fa 10^'* m/s. Assuming vq—IO m/s 
(see, e.g., Ref. p^ ) we estimate the pinning energy as 
Uc = (50 - 55) kT w 0.(^meV for T = 15 K which is 
a quite reasonable valuet3. The other free parameter, 
the sheet current density Jcp, was chosen to provide the 
best agreement of results of flux creep simulations with 
experimental profiles shown in Fig. ^. It is interesting to 
compare the chosen value Jcp = 1.45 x 10^ A/cm with 
the critical sheet current density Jc determined by fit- 
ting the same experimental data to the CSM. We found 
that Jc ~ 0.37 Jcp. For this current density, the effec- 
tive barrier UcS = Uc [{Jcp/ JcT ~ 1] appears substan- 
tially lower than Uc, about Q.lbUc ~ 8kT. It is worth 
noting that variation of /i in the range 1/7-1/3 has mi- 
nor effect both on the calculated profiles B{x) and on the 
relation UeS ~ 8 kT. 



FIG. 8. On the determination of the flux creep parameters 
from the experimental MO data on magnetic relaxation. B{t) 
is the remanent field in the central region of the bridge after 
switching off the external fleld 100 mT. The experimental de- 
pendence ln(— d \nB/dt) versus B'^ with /i — 1/7 is fitted by 
a straight line according to Eq. 

In the relaxation experiment the sample was cooled 
down to 15 K in external magnetic field 100 mT, which 
is about ABc for our bridge. This high field ensures t 
the remanent state is fully penetrated by the current 
Then the field was switched off and the time dependence 
of the remanent field was measured by MO technique in 
the time window 10^-10^ s. From now on we focus on the 
peak field value B{t) observed in the central region of the 
bridge and averaged along its length. Instead of present- 
ing of the relaxation of B directly, we chose to plot the 
logarithmic derivative ln{d In B/dt) versus S^. The rea- 
son for this type of plot is the following. As was shown 
in Ref. |2^, the pinning energy U{x) remains almost the 
same locally at any stage in a relaxation process. That 
is a feature of a self-organized behavior of magnetic fiux 
which is a consequence of the exponential dependence 
V oc exp[—U{J)/kT]. Since U is almost constant v varies 
also slowly in space, and therefore 



dB 
'dt 



dB 
dx 



We can roughly estimate dB /dx as B/w. Moreover, in a 
fully penetrated state dB/dx oc J, and therefore, J (x B. 
This approximation together with Eq. (0) yields 



In - 



dliiB 
dt 



In 



Vo ( Uc 

w \ ki 



const 



B. Results 



Let us first discuss general results of flux creep simu- 
lations. Time evolution of the profiles of current density 
and magnetic field under applying and switching off a 
transport current with density J — 0.26 Jcp is shown in 
Fig. pi The time dependence of the transport current 
was chosen as shown in Fig. |l] in accordance with the 
experimental procedure. For simplicity, only the case of 
zero external field is considered. Since all the distribu- 
tions are symmetric in this case, we discuss the profiles 
for one half of the bridge. 



(a,b) correspond to dif- 
The profiles (1-4) cor- 



The various curves in Figs. 3( 
ferent times as marked in Fig. 1. 
responding to current-carrying states are very similar to 
the ones expected from the CSM. It is interesting to note 
that though the transport current did not change from 
the instant 1 to 2, the curves 1 and 2 in Fig. ^ differ 
substantially. Though both can be described by CSM- 
profiles, they correspond to different values of Jc- In 
particular, the transition from the curve 1 to 2 corre- 
sponds to Jc decreasing in time. Time dependence of 
effective Jc can be seen also from the curves 3 and 4 
for the remanent state after current. A possibility to 
interpret experimental time evolution of magnetic fiux 
profiles by the CSM with time-dependent Jc has been 
previously demonstrated in Ref. 31, An additional fea- 



(8) 



ture of calculated profiles in the remanent state is a peak 
of "negative" current density located not far from the 
bridge edge (cf. with Ref. ^3|) - see curves 3 and 4 on 
Fig. The peak's position corresponds to the annihi- 
lation zone where B{x) changes its sign. This peak is a 
direct consequence of continuity 
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FIG. 9. Results of flux creep simulation for J/ Jcp ~ 0.26 A. Evolution of profiles of the magnetic fleld (a) and the current 
density (b) in the strip plane for the time variation of the transport current shown in Fig. |l| Profiles (a) and (b) correspond 
to different times marked as different points in Fig. |l| 



of the flux flow. Indeed, in the region of low density of 
flux lines their velocity (determined by the ratio J / Jcp) 
should be large to keep the flux flow continuous. 

The dashed line in Fig. |l^ is a result of our simulations 
of the magnetic field profile for / = 0.75 /c at i = 35 ms 
after switching on the current. This time delay corre- 
sponds to the center of the plateau in the time depen- 
dence of the transport current, see Fig. |l|, so it was the 
time when the MO images were recorded. The profile was 
calculated from the simulated current distribution using 
Eq. (^. In a similar way the flux profile for the rema- 
nent state [t = 30 s) was calculated. The result is shown 
Fig. 1^. It is clear that flux creep provides deeper pen- 
etration of the flux into the inner regions in agreement 
with the experiment. 

For the current-carrying state. Fig. the agreement 
is fairly good except two minor discrepancies. First, the 
experimental peaks at the bridge edges are less than the 
ones predicted both by the CSM and by our simulations. 
This is a rather general feature (cf. with Ref. ||) which is 
probably originated from the finite thickness of the indi- 
cator film (in our case 5 /.tm), as well as from imperfection 
of the edges. Another discrepancy observed outside the 
bridge, is obviously related to contact currents. Indeed, 
Eq. (||) is based_ppon the assumption that the exter- 
nal field is smallEil and homogeneous. However, the field 
generated by the current in the contacts is actually in- 
homogeneous in X direction, and far from the bridge it is 
significantly different from the one in its central part. 

For the remanent states after current /, the account 
of flux creep provides good agreement with the experi- 
ments for / < 4.16 A, Fig. |[ However, at larger currents 
the experimental behavior of the trapped flux is qualita- 
tively different. Indeed, both critical state and flux creep 



models predict monotonous dependence of the trapped 
flux on the current. Non- monotonous experimental de- 
pendence can serve as an indirect indication that some 
non-equilibrium process is responsible for the flux distri- 
butions observed after large currents. 
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FIG. 10. Proflles of the absolute value of perpendicular 
magnetic fleld induced by transport current / = 0.75/c at the 
distance h — 8 ^m above the sample. Dots - experiment, 
solid line - calculations along CSM, Eqs. (^),(|l]), dashed line 
- results of flux creep simulation. It is clearly seen that the 
flux creep leads to a deeper penetration of flux comparing to 
the CSM, in agreement with experiment. 
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We believe that formation of the remanent state af- 
ter very large currents is strongly influenced by heating 
effects. As well-known, energy dissipation due to vor- 
tex motion facilitates more intensive motion. As a re- 
sult, mjacrioscopic avalanche-like flux redistributions (flux 
jumps )E3^E3 can take place. Unfortunately, the works we 
are aware of are focused on the flux jumps in the case 
of applied magnetic field and slab geometry. We believe 
that the case of a remanent state in a strip after transport 
current requires a special theoretical treatment. Further- 
more, the flux motion obviously takes place through nar- 
row channels of weak pinning. Consequently, the energy 
dissipation is substantially inhomogeneous that should 
be taken into account in the estimates of local tempera- 
ture. We plan special experiments, as well as the proper 
theoretical analysis as a subject of a future work. 

To get a hint why the heating in the remanent state 
might be different from the one in the current-carrying 
state, we have compared the power dissipation. Accord- 
ing to the estimates, the power dissipation due to vortex 
motion is larger in the current-carrying state. However, 
in the remanent state there is an extra source for dis- 
sipation - vortex-antivortex annihilation. Due to this 
contribution, the total power dissipation in the annihila- 
tion zone can be larger than that in the current-carrying 
state. Consequently the remanent state can be more un- 
stable with respect to temperature fluctuations that the 
current-carrying one. It should be noted that a special 
behavior of vortices in the vicinity of the annihilation 
zone in the remanent state have already been addressed 
in literature. The most unusual feature is meandering 
instability of flux front and its titcbulent relaxation ob- 
served in very clean single crystalscZl. The heat release in 
the annihilation zone is suggested as a probable explana- 
tion of this phenomenon. Another explanation is based 
upon the concept of magnetic field concentration inside 
the bend of a current lineEj. Furthermore, as discussed 
in Ref. |2^, near the annihilation zone, a special behavior 
of flux creep can be expected. It is argued that under 
certain assumptions about _B-dependence of the pinning 
energy, the presence of the annihilation zone destroys the 
normal course of flux creep in the whole sample. 

Another possible source of an instability might be an 
additional heat release in the contact pads. However, ac- 
cording to the estimate given in Sec. HA, the heat release 
in contact regions is negligible. This conclusion is con- 
firmed by the fact that the bridge burn-out took place in 
its central part.E3 



V. CONCLUSION 

Measurements of magnetic fiux distribution in a HTSC 
strip with transport current, as well as in the remanent 
state, are performed by magneto-optical method. The 
experimental results are compared with predictions of 
the critical state model for a strip geometry and with 



computer simulation of flux penetration in the flux creep 
regime. In the current-carrying state, the agreement was 
satisfactory. Simulation of flux creep predicts slightly 
deeper flux penetration than the CSM in agreement with 
experiment. In the remanent state after transport cur- 
rent, the CSM fails to explain the experimental results. 
Our simulations of flux creep allow to describe the flux 
profiles after relatively small current, / < 0.7/c. At larger 
currents the total trapped flux appears substantially less 
than predicted by the flux creep simulations, and it de- 
creases with increasing I. Excessive power dissipation 
in the annihilation zone can be an explanation of these 
experimental results. 

Additional experiments and elaborated theoretical 
models for the remanent state are under development. 
We believe that an important information can be ob- 
tained from time-resolved studies of remanent state nu- 
cleation. The conditions for nucleation of macroscopic 
fiux jumps is also a subject for future theoretical investi- 
gation. 
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